Experimental distillation of squeezing from non-Gaussian quantum states 
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We show theoretically and experimentally that single copy distillation of squeezing from continu- 
ous variable non-Gaussian states is possible using linear optics and conditional homodyne detection. 
A specific non-Gaussian noise source, corresponding to a random linear displacement, is investigated. 
Conditioning the signal on a tap measurement, we observe probabilistic recovery of squeezing. 



^ . 

o ■ 
o ■ 
(N : 

X) ■ 

a^ : 

> 
l> 

o ■ 
(N ; 
o 
^ 
o 

■ 

Oh! 

> ' 

a ' 

a • 

^ : 



X 



PACS numbers: 03.67.-a, 42.50.Dv 

Non-classical states such as continuous variable (CV) 
entangled and squeezed states serve as enabling resources 
for many CV quantum information protocols ll| as well as 
for highly sensitive measurements beyond the shot noise 
limit The efficiency of these appUcations reUes cru- 
cially on the state's nonclassicality (i.e. the degree of 
single- or two-mode squeezing). Therefore, uncontrolled 
and unavoidable interaction of the system with the en- 
vironment and the resultant loss of squeezing in gener- 
ation or transmission should be combated. This can be 
done by using a distillation protocol which probabilisti- 
cally selects out squeezed states from a mixture, hereby 
increasing the output state's squeezing. 

Various protocols exploiting non-Gaussian operations 
to probabilistically distill two-mode squeezed Gaussian 
states have been proposed These protocols are, how- 
ever, experimentally challenging. In addition it has been 
proven that the distillation of two-mode squeezed Gaus- 
sian states by means of more feasible local Gaussian op- 
erations is impossible Q. Similarly, following a set of 
simple arguments, we conjecture that the single copy dis- 
tillation of single-mode Gaussian squeezed states is im- 
possible using only linear optics and homodyne detection. 

There has however been no work devoted to the distil- 
lation of CV Gaussian states corrupted by non-Gaussian 
noise. This occurs naturally in channels with fluctu- 
ating properties, i.e. gain or phase, examples of which 
are the fading channel [5j or channels producing mixture 
noise 0. Recently, non-Gaussian telegraph noise has 
been discussed in qubit systems 7]. Therefore, extending 
the work on Gaussian noise, we pose the question: Is it 
possible to distill single-mode Gaussian squeezed states 
with superimposed non-Gaussian noise using linear op- 
tics and homodyne detectors? We answer this question 
in the affirmative and provide an experimental demon- 
stration. 

The set of non-Gaussian noise sources is large. Thus we 
restrict our attention to a specific case: a squeezed vac- 
uum state perturbed by phase kicks or jitter, attributable 
to either imperfect generation or transmission through a 
noisy channel. Assuming these perturbations cause a lin- 
ear phase space displacement, a convex mixture of two 



Gaussian squeezed states is created 

W{x,p) = (1 - i)Wo[x,p) + iWi{x,p), (1) 

where 7 is the probability of displacement, and the indi- 
vidual constituents of the mixture (i = 0, 1) are described 
by the Wigner functions 
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Here x and p are the amplitude and phase quadratures. 
I^^Xsq and A^Psg are the corresponding variances of the 
input state. X(),p — and xi,pi are the mean values 
of the initial and displaced squeezed states respectively. 
We assume the two individual Gaussians to be equally 
squeezed in x: A'^Xgq < 1, where A'^XgqA'^Psq > 1. 
The first two moments of the amplitude quadrature of 



W{x,p) are 



7x1 



and 



and 



thus the variance of the amplitude quadrature of the cor- 
rupted state of Eq.[T]is A'^X = A^X^g -^7(1 -7)*?. The 
second term here originates from the noise and degrades 
the squeezing. The aim is to recover the squeezing by 
distilling the initial squeezed state from this mixture. 

A schematic of the distillation protocol is shown in 
Fig. n A polarization squeezed state, mathematically 
equivalent to a squeezed vacuum state ^3, ^| , is mod- 
ulated to generate a noisy non-Gaussian state. This 
is incident upon a beam splitter with reflection R and 
transmission T, producing correlated output states. Us- 
ing a Stokes, or equivalently homodyne, detector a given 
quadrature of the tap beam is measured. Conditioned on 
this measurement, the signal is selected only if the out- 
come lies above a given threshold value, as in Ref. 01 ■ 
Due to the correlations between the signal and tap 
beams, the scheme accomplishes a probabilistic distil- 
lation of the noisy input state. A similar strategy was 
proposed to purify decohered Schrodinger cat states . 

We now present a theoretical description of the dis- 
tillation. The tap beam splitter transforms the quadra- 
tures of the Wigner function to, for the transmitted signal 
beam, Xg = VTx + \/Rxu and ps = VTp + VRpv where 
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Xv and Pv are uncorrelated vacuum contributions. We 
write the signal Wigner function after i) detection of the 
tapped signal and ii) post selection of the signal as 



1 



jGi{Xs,Ps)Wi{Xs,Ps)] , 



(2) 



where IT is the success probability and Gi{xs,Ps) is a 
filter function which incorporates the effect of the tap 
measurement and post selection. It thus depends on the 
measured quadrature and the threshold value. Since the 
goal of the distillation is to recover the initial squeezing, 
we consider only the marginal quadrature distribution as- 
sociated with the squeezed quadrature, x. Measuring the 
phase quadrature in the tapped signal pt ^3 j the result- 
ing probability distribution of the squeezed quadrature 
in the signal reads 

P{xs) = ^ {(1 - -/)9oPo{xs) + igiPiixs)} , (3) 

where 11 = (1 — 7)50 + 75i and the individual marginals 
Pq{xs) and Pi{xs) are Gaussian functions with variance 



sq 



RA'^Xy and centered at a;o = and 



Txx, respectively. The filter function in Eq. |31is 

Pth - PiVR 



g, = -Erfc 
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Here pth is the post selection threshold and A^Pt — 
RA'^Psq + TA^Py. After distillation the first two mo- 



ments of the signal are 



{VTxi)/{1 + r) and 



= A^X, + {Txl)/{1 + r) where r = (1 - 7)50/751 
Thus the distilled squeezing is given by 



(l + r)2 



(4) 



The signal variance can be decreased or even the squeez- 
ing recovered by minimizing the second term. The prob- 
ability 7 and the displacement Xi are parameters of the 
noisy process and thus can not be altered in the distil- 
lation optimization. However through the choice of the 
threshold value pth, the ratio between the filter functions 
r can be controlled to yield efhcient distillation, corre- 
sponding to r 00 or r — !■ 0. 

Our experiment for the distillation of corrupted 
squeezed states consists of three parts (Fig. P): the 
preparation, distillation and verification. We observe 
the squeezing of sidebands at 17.5 ± 0.5 MHz relative 
to the optical carrier frequency to avoid low frequency 
technical noise. The preparation of the mixed state of 
Eq. ^ is accomplished by combining a squeezer with a 
controllable noise source. We use a polarization squeezer 
exploiting the Kerr nonlinearity experienced by ultra- 
short laser pulses in optical fibers Using a bire- 
fringent fiber, two quadrature squeezed states can be si- 
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FIG. 1: Schematic of the experimental setup for the gener- 
ation, distillation and verificaiton of non-Gaussian squeezed 
states. 
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FIG. 2: Experimentally measured marginal distributions, cen- 
tered at zero for convenience, outlining the distillation of 
a squeezed state from a non-Gaussian mixture of squeezed 
states; a) tap measurement pt, b) signal measurement Xs. In- 
set: Phase space representation of the mixed state and the 
projection axes used in the measurements. 



multaneously and independently generated. Stable over- 
lap of these pulses allows us to generate Stokes param- 
eter squeezing. Considering the Stokes plane orthogo- 
nal to the classical excitation (6*3, circularly polarized), 
the 'dark' plane {Si — 6*2), it is found that the polariza- 
tion squeezing observed in this mode is mathematically 
equivalent to quadrature vacuum squeezing 0, ^| . We 
treat the two synonymously. The classical excitation can 
then be thought of as a perfectly matched local oscillator. 
From this source we observed A'^Xgq = —3.1 ± 0.3 dB 
relative to the quantum noise level. The anti-squeezed 
quadrature contains the large excess phase noise charac- 
teristic of pulse propagation in glass fibers, here A^Psq = 
-f27±0.3 dB. These noise signals are observed using bal- 
anced detector pairs with 85% quantum efficiency. 
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FIG. 3: Experimentally and theoretically distilled squeezing 
(left) and success probability (right) as a function of post 
selection threshold for two displacements. The threshold is 
given relative the center of the marginal distributions. 



The non-Gaussian noise source is implemented by ex- 
ecuting a fixed phase space displacement of the squeezed 
state with a probability 7 — 0.5. The displacement is 
generated by a phase modulation in one of the linear 
polarization modes at the fiber input using an electro- 
optic modulator at 17.5 MHz. This produces a corre- 
sponding displacement along the S2 polarization after 
the fiber (Fig. inset). The modulation depth governs 
the amount of phase space displacement. By periodically 
switching the modulation on and off, the displacement is 
toggled from maximum to zero at a frequency of 500 kHz. 

This non-Gaussian state, with A'^X = +1.4 ± 0.3 dB, 
is fed into the distiller. It consists of two operations: i) 
the tap measurement of a certain quadrature on a small 
portion of the beam; ii) the signal post selection condi- 
tioned on the tap measurement. The latter could be im- 
plemented electro-optically to probabilistically generate 
a freely propagating distilled signal state. To avoid such 
complications our conditioning is instead based on data 
post selection using a verification measurement. The tap 
and the signal are recorded simultaneously, yielding data 
pairs, and the signal is selected dependent on the tap 
value. These measurements are implemented as Stokes 
measurements in the 'dark' plane. For a circularly polar- 
ized beam the rotation of a half-wave plate introduces a 
relative phase shift between the right-hand circular po- 
larization (squeezed state) and the left-hand polarization 
(local oscillator) when observing the difference signal af- 
ter a polarization beam splitter JdJ. This is equivalent 
to the phase scan of a homodyne detector measuring ar- 
bitrary quadratures in the 'dark' plane. 

The RF photocurrents of the photodetectors are mixed 
with an electronic local oscillator at 17.5 MHz and dig- 
itized with a fast AD converter at 10^ samples per sec- 
ond with a 16 bit resolution. We define our state of the 



electro-magnetic field to be a time window of 1 /Lts. By 
digital filtering and averaging over time bins of 1 fis we 
derive a photocurrent value for each bin. In this pro- 
cess the 1 fj,s time bins of our signal are synchronized 
to the modulator switching period, such that each bin 
is recorded entirely during an 'on' or an 'off' period. 
Thus by measuring the anti-squeezed quadrature in the 
tap on an ensemble of identically prepared noisy states 
we construct the distributions in Fig. [S^a). The simul- 
taneous measurement of the signal beam recorded the 
orthogonal, squeezed quadrature. The modulation was 
chosen such that the variance of the noisy signal was 
just greater than that of the shot noise (Fig.lJfb)). Per- 
forming post selection on this data by conditioning it 
on the tap measurement, we observe a recovery of the 
squeezing. That is, the distilled signal distribution is 
narrower than that of the shot noise (Fig. I2b)). We 
measured A'^Pt = +17.5.0 ± 0.3 dB relative to the shot 
noise. Conditioning on the tap, the noisy signal variance, 
+1.1 ± 0.3 dB, fell to -2.6 ± 0.3 dB after distillation. 

Using the data shown in Fig. 13 the distilled signal vari- 
ance was investigated as a function of the post selection 
threshold. In Fig. El we notice that an increasing thresh- 
old decreases the signal variance, ultimately approaching 
the input squeezing. This agrees well with the expo- 
nential increase in squeezing predicted by Eq. 01 given 
by the dashed line. As the threshold increases, the suc- 
cess probability or amount of distilled data decreases to 
zero causing an increase in the statistical error on the 
variance. Thus a compromise between the post selected 
variance and probability of success must be made. 

The effectiveness of a given threshold depends on i) 
the projection of the displacement onto the measured 
quadrature and ii) the variance of the measured quadra- 
ture. Fig.^Ja) and (b), each with a different threshold, 
shows this effect. The measured tap quadrature was ro- 
tated by an angle /3 (see Fig.|21l, effectively changing the 
displacement size. We observe the best distillation for 
small angles where the displacement {xq — xi) to thresh- 
old {pt) difference is largest. It is seen that the quality of 
the distillation decreases with increasing (3 as the projec- 
tion of the displacement onto the measured quadrature 
increases. We note however that for large thresholds the 
distillation quality is approximately independent of f3 or 
the measured quadrature. 

We also measured the Wigner function, for the first 
time in fiber-based systems, of both the mixed and the 
distilled states. Rotating the half-wave plate in the veri- 
fier (Fig.P) allows observation of all the squeezed state's 
quadratures for a constant tap measurement. We made 
128 equally spaced projections in phase space, each of 
3.5 • 10^ data points to which we applied the inverse 
Radon transformation to derive the Wigner func- 
tions. Fig. EJa) shows the density plots of the Wigner 
function associated with the noisy state; its non-Gaussian 
nature is evident. In Fig. Etb) we present the Wigner 
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FIG. 4: Distilled variance (left axis) and success probability 
(right axis) as a function of the quadrature angle relative to 
the squeezed quadrature (angle /3 in Fig. I^J in the tap mea- 
surement for a post selection threshold of a) 1.3 and of b) 5.3 
shot noise units. 

a) Noisy signal state 



b) Distilled signal state 



FIG. 5: Density plots of the Wigner function distributions for 
a) the non-Gaussian mixed state measured at the verification 
setup and b) the Wigner function of the corresponding dis- 
tilled state. Note: To aid visualization the plots have been 
vertically rescaled by a factor of two. 



function of the post selected signal and observe a dis- 
tribution very similar to that of a single squeezed state. 
Thus the purity of the state relative to the non-Gaussian 
state is increased and a corresponding recovery of the 
squeezing is seen. The shift to the right reflects the post 
selection process as well as the renormalization of the dis- 
tilled data. Further investigation of the Wigner functions 
of fiber-based squeezed states is described elsevifhere |0 . 

We have sucessfully experimentally demonstrated the 
probabilistic distillation of continuous variable non- 
classical states from a non-Gaussian mixture of squeezed 
states. This was accomplished by the thorough investi- 
gation of a specific source of non-Gaussian noise, namely 
a linear shift in phase space. The methods presented 



here can be implemented for many other forms of non- 
Gaussian noise, e.g. phase space rotations, which will be 
the subject of further experiments. Another extension 
of this work presented would be to perform a conditional 
optical operation, i.e. phase shift or displacement, on the 
signal beam. Thus not only the distillation demonstrated 
here but also a purification of the excess noise of the ini- 
tially squeezed states could be implemented, generating 
an even purer non-classical resource than produced here. 
Whilst we have focused on single-mode squeezed states, 
these techniques can assuredly be extended to two-mode 
squeezed systems. This means that continuous variable 
entanglement distillation is possible using local Gaussian 
operations and classical communication if the two-mode 
squeezing is corrupted by non-Gaussian noise. 
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